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ABSTRACT 





Assigning probabilities to RF power picked up on aerospace cables at 
microwave frequencies by using analytical methods was the objective of this 
thesis. The coupling from radiated fields to typical unshielded wires was 
calculated in the frequency range 1 GHz to 10 GHz by utilizing a thin wire 
antenna computer program. The data obtained from the calculations were 
Reed to get the cumulative distribution function of absorbed power and also 
to observe its variation with changing frequency and loading of the wires. 

The results of this analysis compared favorably with those of a previous 
experimental approach. This demonstrates that analytical methods can be 
utilized for calculating and predicting the coupling from radiated fields 


to unshielded wires. 
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me NEED FOR THE STUDY 

In designing aerospace electronic systems to operate in high level 
radiated environments, the susceptibility of electro-explosive devices, 
transistors and integrated circuits to microwave power picked up by con- 
necting cables has long been recognized as a serious problem. Investigations 
into the susceptibility of semiconductor components have shown that they 
respond to microwave energy at power levels of less than a milliwatt. With 
ever increasing radiated environment levels the pickup phenomenon and the 
protection afforded by cable shields have become more important. 

In a radiated environment, circuit wiring can frequently be as efficient 
as tumeddipole antennas. This provides much more than a few milliwatts of 
power to cable terminations when wiring is exposed to high radiated 
environments. 

Designing aerospace electronic systems to survive and operate in high 
power radiated environments requires techniques for predicting and evalu- 
ating the coupling from radiated fields to susceptible components. This 


permits system designers to establish the overall required system protection. 


B. STATEMENT OF THE PROBLEM 

Inefficiently shielded aerospace cables may pick up enough energy to 
cause the electro-explosive devices of a missile to prefire or the missile 
to change itS course. Even with good shielding, pickup sufficient for deto- 
nation can occur in a severe environment. 

Therefore in missile design, for successful operation in hostile envi- 


ronments, protection of electronic circuits against high power radiated 


fields must be provided. 
10 





The necessity of evaluation and prediction of RF coupling to permit 
system designers to establish the overall required system protection was the 


m@eivine force for this thesis. 


C. PREVIOUS WORK 

In an electronic circuitry, a cable can be thought of as an antenna which 
couples power from the radiated field to the sensitive component and an ef- 
fective aperture can be used to describe the coupling problem. 

McDonnell Douglas Astronautics Company - East (MDAC-E) has worked on 
solving the high frequency cable coupling problem, and published the results 
of its studies as a menor eS According to this report the studies have 
resulted in an approach which divides the problem into two parts: (1) ef- 
fective aperture of the unshielded cable and (2) shielding effectiveness of 
the cable shield. These can be stated in equation form: 


P.(W)=P, (W/m?) x A, (m*) x SE 


Bere: P. = power received by cable termination 
tgum power density of radiated field 
A, = effective aperture of unshielded cable 
SE = shielding effectiveness of cable shield 


MIL-STD-1377 (NAVY) provides a good method to evaluate shielding para- 
meters independent of other variables (polarization, aspect angle, cable 
length, etc.). Due to the complexity of practical shielding and enclosure 
discontinuities, shielding effectiveness is determined by empircal methods 
in this approach. However, for determining the effective aperture in MDAC-E 
approach both analytical and experimental methods may be utilized. 

To validate both the MDAC-E approach and the MIL-STD-1377 shielding 


effectiveness measurement technique, data have been processed from a large 


il 





number of anechoic chamber measurements of RF pickup at 3.0 GHz and 9.1 GHz 
for shielded and unshielded twisted wire pairs. The cumulative distribution 
of RF pickup for each test configuration has been calculated. The data were 
put in a form suitable for assigning probabilities. Results of this work 
showed that the MDAC-E approach is valid and that MIL-STD-1377 provides an 
accurate method to measure shielding effectiveness. The data also indicated 
that the maximum effective aperture of an electrically long unshielded twisted 
wire pair can be approximated by a half-wave dipole at microwave frequencies. 
It was found that the unshielded cables had a measured maximum effective 
aperture very close to that calculated using a half-wave dipole model. 

In statistical analysis of MDAC-E approach an antenna pattern of received 
power was plotted while the cable was exposed to a plane wave in an anechoic 
chamber and rotated 360 degrees in azimuth. A separate pattern was recorded 
every 5 degrees in evaluation, from 0 through 180 degrees. 

Using these data the cumulative distribution of RF pickup was obtained. 
Figure I-l shows the cumulative distribution of the data points for 9.1 GHz 
and the 360 degree ealeiaa cable and the unshielded cable which had the same 
terminating impedances. The two curves are very nearly equidistant to all 
points. If one could measure or calculate the coupling to an unshielded 
cable, then shield termination techniques can be evaluated using MIL-STD-137/7, 
and RF coupling to the shielded cable can be predicted. 

However a point-to-point corresponding reduction in RF coupling does 
not exist. That is, the antenna patterns for the shielded and unshielded 
case are not similar in a geometrical sense. In actual use, the cable will 
be placed in proximity to objects which will change the RF coupling for a 


particular aspect angle considerably, but the statistical distribution will 
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cable. Field intensity=1 V/m 
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not change. Therefore, RF coupling to the shielded cable must be predicted 


Sieaestatistical basis. 


DeeeoCOPE OF THE STUDY 

In this thesis study the cumulative distribution of RF pickup of unshield- 
ed wires was evaluated by utilizing completely analytical methods. Since 
unshielded wires can be thought of as antennas, antenna modeling programs 
can be used to calculate the RF power pickup by computers. 

M B Associates Antenna Modeling Program (AMP) was used to calculate the 
currents on a straight wire in a free space radiated environment. The mag- 
nitudes of currents for different Rect and polarization angles provided 
the data to compute the cumulative distribution of RF pickup energy. For 
the purpose of computing the RF power pickup and plotting the cumulative 


distribution another computer program was used. 
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It. ANALYSIS AND PROGRAMS 


A. MODEL OF STRUCTURE AND STATISTICAL ANALYSIS 

The test structure used by MDAC-E was a model of an integrated circuit 
NAND gate driving another NAND gate through a twisted wire pair exposed to 
a microwave field. The IC impedances were measured and simulated using 
impedance matching networks so that a receiver with large dynamic range 
could be used to monitor power received by the termination. The test cables 
were a pair of twisted wires, one twist per inch, 23-5/8 inches long. One 
conductor was connected to the center pin of type N to type BNC adapter at 
each end. The other conductor was soldered to the connector shells. With 
this test structure several sets of measurements were made for different 


configurations, three of these configurations are shown in Table II-l. 


CONFIGU- FREQUENCY LOADS 

RATION TERMINATION CALCULATION AT 
1 oF imGliz ge fs Oi 30 
2 3.0 GHz 40-523 24-428 
3 iO =GHz 50 50 


Taplesti-l;  Ceantiguuacion record 


In this study, the test structure explained above was simulated to cal- 
culate RF pickup power and its cumulative distribution calculated by computer 


programs. In order to be able to compare these results with those of MDAC-E, 
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all the parameters of the model of the test structure were chosen to be the 
same as MDAC-E's parameters shown in Table II-l. 

imestest structure which consisted of a pair of twisted wires and tWo 
loads connected to their ends was modeled by a single wire having one load 
at each end. Since one of the wires in the test structure was used as ground 
cable, the current produced by incident plane waves on this cable flowed to 
" ground and made no contribution to the amount of power picked up by the ter- 
minating impedance. Therefore the ground cable of the test structure was 
omitted in the model; 

The 23~5/8 inch length of the model wire was kept constant and several 
positions of the loads were tried for optimum correlation with MDAC-E's result. 
It was found that the best location for the loads was the ends of the wire, 
as it was in the real case. These positions of the loads gave the optimum 
correlation for both configuration 1 and configuration 2. The comparison 
for configuration 3 was not done, since the cumulative distribution of RF 
pickup power obtained by MDAC-E was not available. 

Varying the load positions demonstrated that the RF pickup power was 
strongly dependent on the load positions of the structure. The effect of 
load position on the cumulative distribution of RF pickup power is shown in 
Figure II-1 for configuration 2. The plotted distribution curves are for 
positions at 0.05 wave length and 0.15 wavelength from the ends. Variation 
of load locations of 0.1 wavelength caused the cumulative distribution <t& 
shift about 4 dBm. However further study showed that the amount of shifting 
was not directly proportional to the variation in load positions but it 
could fluctuate randomly between 10 and 2 dBm. 

The model of the test structure in three-dimensional coordinate system 


is shown in Figure II-2. Load 25 was used as termination and power and cumu- 


lative distribution calculations were made for At 
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In the statistical analysis, the first step was calculation of currents 
at load zy for incoming plane waves which had different aspect and polari- 


zation angles()).This was done by uti’ ’»iug MBA's Antenna Modeling Program. 


A 





XY Figure II-2 
Model of the test structure in three-dimensional coordinate 


system. 


The elevation incidence angle and polarization angle were varied in 
5-degree increments in calculating the spherical receiving pattern of the 
Structure. Since the structure was symmetrical in POH Ome CASR the 
aximuth angle (4) of incident plane wave was kept constant. 

To determine symmetry of polarization angles, currents were evaluated 
by varying both theta and eta angles from 0 to 180 degrees. It was found 
that incident plane waves having polarization angles between O and 90 degrees 
produced the same currents in magnitude as the plane waves having polariza- 
tion angles between 90 and 180 degrees, but their phases were 180 degrees 
different. Therefore the eta angle scanned only 9 to 90 degrees. Since 
the model of structure was unsymmetrical in see direction, a 180 degree 


sector of 9 angle was covered. 
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B. CALCULATIONS 

Calculations made to obtain the cumulative distribution of RF pickup 
power can be analyzed in two groups: (1) evaluation of currents at load zy 
produced by incoming plane waves and (2) evaluation of power in dBm, calcu- 
lation and plotting of the cumulative distribution. 

Two separate computer programs were utilized to achieve the two groups 
of calculations stated above. 

Evaluation of currents at load Z, was done by the AMP program. It is 
based on the thin wire electric field integral equation which relates the 
exciting electric field to the induced currents on some specified thin wire 
geometry. The integral equation is reduced to an N-dimensional system of 
linear equations by representing the current in terms of N sinusoidal basis 
functions and enforcing the integral equation at N discrete points. The 
unknown current vector is related to the excitation vector by a coefficient 
matrix which essentially depends only on the geometry of the structure. The 
unknown currents in the matrix equation are solved for by using the Gaus- 
Doolittle elimination procedure. Fields patterns can be calculated by sum- 
ming the contributions of the N currents. The basic elements of the program 
used for solving the electromagnetic characteristics of a structure from the 
electric field integral equation are relatively simple. However, the exis- 
tence of impedance loading, ground effects, and network options in the progran 
increases its complexity. 

Another computer program was written to calculate the RF pickup power in 


dBm and its cumulative distribution. It utilizes subroutine DRAWP to plot 


the cumulative distribution of RF pickup power. 
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Since it was impossible to indicate the zero power point on the dBm 
scale, all zero magnitude currents were extracted from the data. Hence the 
number of data points were reduced from 703 to 630. It was assumed that the 
probability of having each RF pickup power value from 630 possible values 
was equal. In order to calculate cumulative probabilities, the number of 
RF pickup power values which were greater than some certain specified level 
was counted and divided by 630. The Antenna Modeling Program calculates the 
currents on a structure for an incidence electric field of 1 V/m. In MDAC-E 
method, measurements have been made for an electric field of 19.4 V/m. To 
be able to compare the results of this study with MDAC-E's results, current 
magnitudes were multiplied by 19.4 to give the current magnitudes for an 
electric field of 19.4 V/m. The program and its flow chart are shown in 


Appendix A. 
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ITI. RESULTS OF CALCULATIONS 


The cumulative distributions of RF pickup power for the first and second 
configurations in Table II-1 are shown in Figure III-1 and Figure III-2. In 
order to compare the results of MDAC-E's experimental approach and the com- 
putational method employed in this study, they were plotted on the same coor- 
dinate system. In both figures, it can be observed that there is a good 
correlation between experimental and analytical curves around the 50% cumu- 
lative probability point. However the difference between two curves increases 
towards the high and low power ends. This deviation of analytically obtained 
curve from MDAC-E's curve was caused by the differences between the actual 
test structure and the model of it. These differences could be divided into 
two groups: 

1. In the test structure the two loads were connected to the ends of 
twisted wire pairs. But, because of the limitations of Antenna Modeling 
Program it was impossible to place the loads at the end points of the struc- 
ture model. AMP makes its calculations assuming loads are at the center of 
segments. This situation causes a small part of the wire to appear after the 
load at the ends of the structure model. Since the current magnitudes dimi- 
nish towards the end of structure, the measured values of power were smaller 
than the calculated values due to the larger distances of loads of structure 
model from ends. This phenomenon affected high power end of the analytically- 
obtained cumulative distribution and caused it to have larger values than the 
experimentally-obtained curve. 

2. In order to model the test structure, solid wire was used instead of 


Svteeea wire. Since the effective cross-section area of twisted wires is 
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Comparison of experimentally and analytically obtained 


cumulative distribution curves for configuration 1 
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larger than solid at radio frequencies, this also caused a small amount of 
difference between the two distribution curves. 

The cumulative distribution of RF coupling for the third configuration 
is shown in Figure III-3. Comparing the three plots, it can be seen that 
the shape of distribution curves were indepndent of parameters used. How- 
ever there was a shifting in the negative direction and an increase in the 
slope of the curves with increasing frequency. The frequency of the radiated 
environment and the positions of the loads on the structure were the major 
factors which change the amount of RF coupled energy. The values of the 


impedances didn't change the cumulative distribution appreciably. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Assigning probabilities to RF power coupled from radiated fields to 
aerospace cables by utilizing computational techniques was the purpose of 
this thesis. This objective was accomplished by use of two computer prog- 
rams. A test structure which was used previously by MDAC-E to obtain the 
same objective experimentally was modeled to calculate the current magni- 
tudes produced by incident plane waves. Current magnitudes calculated by 
use of MBA's Antenna Modeling Program provided the data for another computer 
program to calculate and plot the cumulative distribution of RF pickup power. 

In modeling, choosing the load positions of the structure model the same 
as the load positions of the structure to be modeled gave the best correla- 
tion with MDAC-E's result. Also a single wire was utilized to model the two 
wire pair of the test structure. 

The results showed that the frequency of the radiated environment and 
the positions of the loads on the structure were the major factors which 
affected the amount of coupled RF power. The effect of the other parameters 
on coupled RF power was relatively small. 

Comparisons of the resulting curves of this study and MDAC-E's experimen- 
tal approach was made and it was found that there was good correlation be- 
tween them. Around the 50% cumulative probability point, agreement between 
the two distribution curves was even better. That is, evaluation and predic- 
tion of RF coupling from radiated fields to aerospace cables can also be done 
by utilizing computational techniques, which are less time consuming. Many 


types of wiring found in aerospace electronic circuitry can be modeled by 
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these computational techniques. Knowing the coupling to unshielded cables, 
the MIL-STD-1377 (NAVY) method of measuring the shielding effectiveness of 
shielded wires can be used to accurately predict RF coupling from a radiated 
field to the terminating impedance of shielded cables. 

In practice, there is more than one subsystem which is suceptible to a 
radiated environment. The probability that at least one of the subsystems 
is vulnerable increases with the complexity of a system. For a more realis- 
tic prediction of RF coupled power to aerospace electronic systems, the 
variation of cumulative distribution of RF pickup power with the number of 
subsystems should be determined. This may be an area of possible improve- 
ment of analytical methods utilized in this study to permit system designers 
to predict the overall system protection. 

Another consideration is to calculate the probabilities of malfunction 
and burnout of the susceptible components of aerospace electronic systems as 
a function of distance from a representative high power transmitter or as a 
function of the power level of a radiated environment. This is looking at 
the same problem from another point of view and may be an alternative compu- 


Pational solution. 
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COMPUTER PROGRAM TO CALCULATE AND PLOT 


CUMULATIVE DISTRIBUTION OF RF PICKUP POWER 


APPENDIX A; 
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Flow chart of the program which calculates and plots the 
cumulative distribuciaun J1 wt pickup power 
OIMENSION CURENT (630) »,POWR( 630), POWRM( 630) »06M(630) 


DIMENSION DELDB(340),PROB(340),URENT(630) 
INTEGER*4 ITB(12)/12*0/ 


REAL*4 RTB(28)/28*0.0/ 

Piet3s), =5 

IT6(4) =7 

RTS(1) =20.0 

RTb(2) =0.2 

RL =50.0 

CKOL =630.0 

KOL =630 


meeawew @Ee am abPpeanw fb  @= 42 OE ee SE ee OS SS ee ae ean a EEew ae ES a a eee Se ee eee eee ee ee 


**¥*WRITE (6911) 


el) FORMAT(C'1"',48X,*POWER IN UNITS OF DBM‘) 
**xxXWRITE (6,13) 
3 FORMAT(C' a 


CURRENT MAGNITUDES ARE READ 


*¥**READ(5S,35) CURENTCI) »f=1,KCOL) 
35 FORMAT(¢(iX,D10.4) ) 
THIS COOP CALCULATES THE POWER IN UNITS OF DBM 


+ DO + 
H#t¢ttt¢ttttt 1 + 
+ [T=1,KOL + 
t—------------ + 
CURENT(I) =19.4C9544*URENT(T) 
POWR( IT) =RL*¥(CURENT(1)**2) 


qn seeaw @] So aw ® ] =e ae eanw Se eweweawens SF S oe SS SS SBS SS Se es SE eT eT Tee Ue eee lc lc ee leer eer lel ee oe 


POWER [IS CONVERTED FROM UNITS OF WATS TU MILLIWATS 


an =e 2 aw 2 2 aga SF ewe—s ae ew es a Ss ee Pee ee ee ee ee eee eee 
—_— ow ae ae om oe ewe ae ee a ce ei se iiss ees ae ae ae ie ee Se oe 


~~ ep ow ae ee eee a ieee ae ee ap ap ee ee cee ame ee ee So eee ee ee ee eee 


HHH HHH EH He HTH He He ete este 


—_ew GE ee SO aw ew aw ew ew aoe ee ewe es ee i ee ee ee eee ie ee 
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POWERS ARE WRITTEN IN UNITS OF DBM 


*XXWRITE(6,56) (DBM(1),T=1,KOL) 
Se et het) 


—_—~ au ca car «= ae = oe cep ee cme cm oe ae eo eee ee ee ee Ss om cm ee a ee ee ee ee es ee ee eee eee ee 


**XWRITE (6712) 

FORMAT ('1*,48X,'CUMULATIVE PROBABILITiES') 
**X*xWRITE (6,13) 
THIS LOOP CALCULATES THE CUMULATIVE PROBABILITIES 


me cee es ee eee ee ee + 

+ DO + 
ttt ttt ett 7 + 
+ + K=1,340 4 
+ $t------------ + 
: | 
+ 
ee eee ee eee eee ee ee ~~ ewes Bw ew wwe es = = 
+ { CONTER =0.0 | 
ra > ame i em ee es es ees ee ee ee ee ee ee ee ee es ee ee ee ee ee ee ee es ee ee ee wee 
: | 
+ 
$$ mee ee een ne ee a 
. + DO + 
Ftt¢ttttttt 4 
+ + T=1,KOL + 
$e ee + 
t 
+ 
+ 
+ Abe: ome 
+ es HE * 6 
Re Cees ELE TK) x 
+ s @ Keen een ee 
+ ~ g . * T ] CONTER=CONTER+1.9 
+ Bg Ke ee ee ee eH 
. 
+ F 
$e a a ae ee ee a ae ee ee ee 
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ii 


35 


SHEEEEEEEHREEHE EEE EEE EEE EHS 





t¢+tttttCONTINulE 
CUMULATIVE PROGARBILITIES ARE WRITTEN 


ESN SS ad (PROB (K) »K=1,340) 


| 
ee ee gel bee) 


oan ee Oe ee ee ee ew ee eo See 2 ew es Sw ew ew 2 2 SO Oe ow ow 2 ow Sw ew es is es ee sae es a 


CUMULATIVE PROBABILITY VS. POWER IS PLOTTED 


CALL DRAWP(340,DEL0B,yPROB;,1TB;RTB) 
STCP 


Sil 





Symbol dictionary: 


external routine 


ALOG10 


magnitude of current at load Z, for electric field of 19.4 V/m 


CURENT 
CKOL = number of current magnitudes (real number) 

DBM = power at Z, in units of dBm 

DELDB = variable name of increasing value of received power 


DRAWP external routine 


ITB = array in subroutine DRAWP containing miscellaneous information for 
the function to be plotted 

KOL = number of current magnitudes (integer number) 

POWR = power at load Z, in units of watts 

POWRM = power at load Z, in units of milliwatts 

RL = resistive component of load Z) 

RTB = array in subroutine DRAWP containing information for scaling, and 


producing a title for the completed graph 


URENT = magnitude of current at load Z, for electric field of 1 V/m 
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APECHNDiwep. EFFECT OF THE MEASUREMENT CABLES ON PICKUP POWER 


In determining coupled RF power from incident plane waves to susceptible 
components of electronic circuits by measurement, the cables between meters 
and susceptible components also pick up radiated power and cause inaccurate 
power readings on the meters. In order to observe the effect of measurement 
cables on RF power absorbed by a susceptible component in an experimental 
method, the structure model shown in Figure B-1 was analyzed with and without 
measurement cables. In modeling and calculations, the methods explained in 
II-A were utilized. Hence the measurement wire pair was modeled as a single 


wire. 





Measurement 
cable 


Figure B-l 


Model of the test structure and its measurement cables 


Coupled power was calculated for ten different azimuth, elevation and 
polarization angles. The results of these calculations and the parameters 
of the structure are shown in Table B-l. It is seen that for some aspect 


and polarization angels the difference in received power for the structure 
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models with and without measurement cables was appreciable. This difference 
was large especially for incident plane waves having polarizations parallel 
to measurement cables. Therefore in MDAC-E's experimental method which was 
utilized to determine cumulative distribution of RF pickup power, the measure- 
ment cables could be a source of error. However, since the calculations were 
made for unshielded measurement cables in this study, the difference due to 


shielded measurement cables must be expected to be much less. 


34 





Frequency L Terminated imnedance 
3.0GHz 23-5/8 inches 24-328 


POWER CALCULATIONS 


4 ) ») Power with con- Power without Difference 
nection cables connection (dBm) 
(dBm) cables (dBm) 

0 30 0 =O. 2740 Oo 77 0.0763 
45 45 0 -11.2959 -13.5604 2.2645 
45 45 4S =f) US -8.9120 Loyal 
ao 6 240 60 = oOZ -6.9389 2475 
ZO 6S e210 30 Sis yisSh: -9,6176 tea LAs) 
nO. 55 10 is SSI: 11.4176 eS 7 9 
70 o3 80 -4,4082 = ls MES OSA 
65 v5 50 -1.8102 las ROS 
20 40 40 -9,.4609 -8.8612 0.5997 
90 0 0 SAO 50.0205 48.4827 


Table B-1: Parameters of structure model and calculation 


record 
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